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FOREWORD 


The report describes the development of the Tridiagonal Reduction Method 
for real eigenvalue analysis and the implementation of this technique in NASTRAN. 
A follow-on report. Part II, will subsequently be issued for a complex eigenvalue 
version of the method. 

All major steps in the analytical development are documented. In addition, 
a detailed summary of the computational procedures, explanatory flow diagrams and 
user data-preparation instructions are provided. 
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l. INTRODUCT ION 

The Tridiagonal Reduction or PEER Method Is an automatic matrix reduction 
scheme whereby the eigensolutions in the neighborhood of a specified point in the 
eigenspectrum can be accurately extracted from a tridiagonal eigenvalue problem 
whose order is much lower than that of the full problem- Specifically, the order, 

m, of the reduced problem is never greater than 


m = 2q + 10 


Cl) 


where q is the desired number of accurately computed eigenvalues. Thus, the 
intrinsic power of the method lies in the fact that the size of the reduced 
eigenvalue problem is of the same order of magnitude as the number of desired 
roots, even though the discretized system model may possess thousands of degrees 
of freedom. The process is effected without arbitrary lumping of masses or 
other physical quantities at selected node points and thus avoids one of the 
basic weaknesses of the Guyan Reduction Method (Reference [1]) and other tech- 
niques (References [2] and [3]) requiring a judicious selection of the degrees 
of freedom to be retained. 

Tridiagonal reduction was first suggested by Crandall (Reference [4]) as 
a truncated version of the Lanczos Algorithm (Reference [5])- However, it was 
soon discovered that the original scheme possessed numerical instabilities 
(References [6] and [7]). The necessary improvements to correct these weak- 
nesses were made by Ojalvo and Newman (Reference [8]) who were the first to 
develop a successful tridiagonal reduction program for large scale structural 
vibration problems. Further refinements were later introduced by Newman and 
Pipano in the PEER computer program (References [9] and [10]), including the 
following extended features: 

1. Highly efficient numerical computation schemes which take 
advantage of matrix banding and sparsity. 

2. Calculation of upper and lower error bounds on the extracted 
eigenvalue estimates. 
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3. 


Accommodation of singular mass matrices and stiffness matrix 
singularities associated with rigid body modes. 

The last capability, which is also present in the Inverse Power Method 
with Shifts, overcomes a basic restriction of the Tridiagonal (Givens) Method, 
namely, the need to eliminate massless degrees of freedom- 

From the standpoint of computational speed, the Tridiagonal Reduction 
Method is almost as fast as the Givens and Householder methods (References [11] 
and [12]) when all the eigensolutions are calculated, and becomes increasingly 
more efficient as the number of required eigensolutions is reduced. In addition, 
in order to avoid prohibitively long running times, both the Givens and House- 
holder methods require the use of a relatively large computer central memory 
for even moderate problem sizes, while the Tridiagonal Reduction Method is ex- 
tremely efficient with regard to core requirements. 

As shown in Section 2.1, the Tridiagonal Reduction Method employs only 
a single initial shift of eigenvalues and hence usually requires only one matrix 
decomposition. It consequently tends to be much more efficient than the Inverse 
Power Method when more than one or two eigensolutions are required. 

This report describes the development of the Tridiagonal Reduction Method 
and its implementation in NASTRAN for real eigenvalue analysis as typified by 
structural vibration and buckling problems. 


The restrictions on the use of the method for real eigenvalue analysis in 
NASTRAN are as follows: 

1. For structural vibration mode applications the method extracts 
a preselected number of eigenvalues which are closest to a 
specified shift value, X , rather than computing the eigenvalues 
in a prescribed range. 

2. In buckling problems, a preselected number of eigenvalues of 
smallest magnitude are obtained, i.e., no shifting is performed. 
Physically, this implies that the buckling load parameters, 
whether positive or negative, are computed in order of increasing 
magnitude. 
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The basic steps employed in the method are as follows. First, the 
initial eigenvalue problem 


[K - X M] {<)>} = 0 (2) 

is converted to a symmetric inverse form 

[B] {X} = A[D] {x} (3) 

where 



and X is a shift value which is used only in structural vibration mode 
o 

applications. Second, the tridiagonal reduction algorithm is employed to 
transform Equation (3) into a tridiagonal form of reduced order. Third, the 
eigenvalues of the reduced matrix are extracted using a Q-R algorithm similar 
to that described in Reference [13]. Fourth, the corresponding eigenvectors 
are computed and converted to physical form. Finally, upper and lower error 
bounds on the extracted eigenvalues are obtained. 

The development of the method is set forth in Sections 2 to 6. A 
detailed summary of the computational procedures used in NASTRAN, explanatory 
flow diagrams and user data preparation instructions are provided in Sections 
7 and 8. In addition, user information and error messages and optional diag- 
nostic output relating to the Tridiagonal Reduction Method are described in 
Section 9. 
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2 . 


THE TRIDIAGONAL REDUCTION METHOD 


2.1 Preliminary Operations 

The problem is to find a specified number of real eigenvalues and corres- 
ponding eigenvectors for 

[K - X = 0 (1) 

a 

It is further required that these eigensolutions constitute the set lying 
closest to a specified point, A q , in the eigenspectrum. 

The definitions of the eigenvalue, A , the matrices [K] and [M] , and 

a 

their mathematical properties, depend on the type of problem being solved within 
the NASTRAN environment. For real analysis, which is the subject of the current 
report, only two separate problem types need be considered; structural vibration 
and buckling problems. The matrix definitions and mathematical distinctions for 
these two cases are summarized on the following page. 
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TABLE 1 - Problem Formulations 


Problem 

Type 

Quantity 

Definition 

NASTRAN 

Notation 

Most General 
Properties 

Structural 

Vibration 

Modes 

[K] 

Stiffness Matrix - 
analysis set 

[K ] 
aa 

Symmetric, non- 
negative, semidef inite 
matrix 

[M] 

Mass Matrix - 
analysis set 

[mi 

aa 

Same 

X 

a 

Square of a 
circular natural 
frequency 

2 

to 

Positive 

Buckling 

[K] 

Stiffness Matrix - 
analysis set 

ik i 

aa 

Symmetric , positive- 
definite matrix 

[M] 

' 

Differential Stiff- 
ness Matrix - 
analysis set 

J 

aa 

Symmetric , indefinite 
matrix 

X 

a 

1 

Buckling Load 
Parameter 

-X 

Positive or negative 

! 


The essential mathematical differences between the two types of problems 
center around the properties of the [M] matrix, which is non-negative for vibra-^ 
tion mode problems* but indefinite for buckling problems, thereby permitting the 
existence of both positive and negative eigenvalues in the latter case. In addi- 
tion, the stiffness matrix may .be singular for vibration problems while it is al- 
ways positive definite in buckling applications, which implies that the buckling 
analysis is performed on a kinematically stable structure. 

In summary, the two problems under consideration are of the forms 


and 


[K - 0) M ]{<{>} » 0 
aa aa 


(structural vibrations) 


(2a) 


[K + AK d ]{<(>} = 0 (buckling) 

aa aa 


(2b) 


Further, if the user requests vibration inodes in the neighborhood of a 
specified frequency, u > q . Equation (2a) can be written as 
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(3) 


[K]{<j>} = X' [M 

where 

[K] = [K - A ] (4a) 

aa o aa 

and 

X’ = w 2 - CJ 2 (4b) 

o 

The resulting effective stiffness matrix, [K] , is indefinite in this case, since 
it possesses both positive and negative eigenvalues. This requires that a non- 
square root decomposition scheme be used in subsequent operations. However, 
o> o = 0 is taken as a default value, or it may be specified by the user. In this 
case, a specified number of natural frequencies starting with the lowest will be 

computed. In order to utilize a more efficient Cholesky decomposition of [K] 

2 

under these conditions, a small negative shift X^ = -a (see Section 5) is used, 
yielding 

[K] = [K + a 2 M ] (5a) 

aa aa 

and 

X 1 = to 2 + a 2 (5b) 


It is easy to prove that the resulting effective stiffness matrix [K] is 
positive definite provided that the system masses generate positive kinetic 
energy due to any kinematically admissible rigid body motions of the structure. 
This requirement is always satisfied by the mass matrix in a physically well 
posed problem, thereby allowing a Cholesky square-root decomposition to be per- 
formed when the roots are computed in the neighborhood of zero. Since no shift- 
ing is performed in buckling problems, the effective stiffness matrix is 
[K] = [K ] , which is always positive definite, again permitting the use of a 
Cholesky decomposition. 


In any event, a decomposition or factoring of [K] is next performed: 


or 


[K] = [L]['dJ[L] T 

(shifted vibration mode problems) 

(6a) 

[K] = [C][C] T 

(buckling problems or vibration 

(6b) 


modes in the neighborhood of 



zero desired) 
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•where [L] and [C] are lower triangular factors and [“'•dj is a diagonal matrix. 

To facilitate computation of eigenvalues closest to the point of interest 
within the eigenspectrum, inverse forms of the eigenvalue problems are employed, 
as in the Inverse Power Method with Shifts. 

The general form of the inverse problem may be written as 

[B] {X} = A[D]{X} (7) 

where the above terms are defined as follows: 


TABLE 2 - Inverse Eigenproblem Definitions 


Problem 

Type 

[B] 

[D] 


A 

1. Shifted 
Vibration 
Modes 

[M aa ][L " 1]T[ ' d ' ] " 1[L " 1][M aa ] 

[M aa ] 

{*} 

1 

2 2 
0) -0) 

o 

2. Unshifted 
Vibration 
Modes (in 
the neighbor- 
hood of zero 
frequency 

[c 1 3 [m ] £ c 1 3 T 

aa 

1 

[I] 

(Identity 

Matrix) 


1 

2 2 
w +a 

3. Buckling 
Modes 

I c_ 1 3 [K d 1 [C _1 ] T 
aa 

[i] 


1 

X 


The above triangular matrix inverses are treated as purely operational 
symbols, since in actual numerical computations vectors defined, for example, by 


{a} = [L 1 ]{b} (8a) 

{a} = [L _1 ] T {b} (8b) 


are obtained from the solutions of 
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II 


[L] {a} = {b} 

(9a) 

[L] T {a}= {b} 

(9b) 


employing forward and backward passes. 

2.2 The Reduction Algorithm 

A reduction of the order of the eigenvalue problem. Equation (7), is 
effected through the transformation 

{x} = [V] {y} (10) 

nxl nxm mxl 

where {X} is an approximation of {x}, n is the order of the unreduced problem, 
and m n. The transformation matrix is taken to be orthonormal to [D]*, so 
that 

[V] T [D][V] = [I] (11) 

From Equations (7), (10), and (11) it is seen that 

[A] {y } = A{y} 

where 

[A] = [V] T [B][V] 
mxm 

and A is an approximation of the eigenvalue, A. 

Thus, Equation (12) is an mth order eigenvalue problem where m <_ n. 

The value of m is established according to the criteria given in Section 2.3. 

The essence of the reduction scheme lies in the choice of the transforma- 
tion matrix [V] . In the present case the Lanczos algorithm is used to build up 
the [V] matrix, vector by vector, i.e., 

*In problem types 2 and 3 (see Table 2), [D] is the identity matrix and, 
therefore, [V]T[V] = [I] . 


( 12 ) 

(13) 
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( 14 ) 


[V] = [{v^, (v 2 > tv^}], 

nxm 

such that the reduced mxm matrix [A] is tridiagonal and its eigenvalues 
accurately approximate the roots of Equation (7) having the largest magnitude 
(or, equivalently, the roots of the physical model closest to the specified 
point of interest in the. eigenspectrum) . 

Define the matrix 


[B] = [D hfB]*, 


(15) 


and let 


d i+i {v i+i } = { Vhi } = tS]{v i } “ a i,i {v i } ‘ a i,i-i {v i-i } 


- a . .{v,} 
x,l 1 


= [B]{v } - l a (v }; 1=1, m-1 

j=l ,J 3 


(16) 


where {v^} is a starting vector (see Section 3) and all m {v}-vectors are 
orthonormal to [D] . This implies that 


a. . = {v.} [B] (v . } 
1.3 3 a. 


(17) 


while d^ + ^ is a normalizing factor given by, 


d i+l = [{ ^i+l }T[D]{v i+l }]1/2 


(18) 


The recurrence relationship. Equation (16), when carried out to one addi- 
tional step (i.e., over the range i = 1, m), can be expressed in the following 
matrix form: 


*For shifted vibration mode problems, where [D] = [M aa ], this operation is. purely 
symbolic (since [M aa ] may be singular) and Equation (15) merely implies that 

m = i. 


aa 
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( 19 ) 


[5) [V] - [V] [HI + d^IG^] 

nxn nxm nxra mxm nxm 

where [H] is an upper Hessenberg matrix given by. 



a il 

a 21 

a 31 

• » • • a - 
ml 


d 2 

a 22 

a 32 

• * • * a m2 

[H] = 


d 3 

a 33 

• 

• 





• 

* * 

d a 

m mm 


«W ■ l{0} (0) ! : <°> : (Vi 11 - 

T 

Premultiplying Equation (19) by [V] [D] , it can be seen that 

m T [B][V] = [V] T [D][V][H] + [0] 
or 

[A] = [H]. 


( 20 ) 


( 21 ) 


( 22 ) 

(23) 


However, the matrix [A] is symmetric and hence [H] is symmetric, requiring that 
it be tridiagonal. It follows that 


and 


a. . = 0 ; j < i-1 


d. = a. . , 

l l, i-l * 


while Equation (19) assumes the form 


(24) 

(25) 


[B][V] = [V] [A] + d m+1 

As a consequence of Equations (23) to (25) 
value problem. Equation (12), is 


■ < 26 > 

, the reduced tridiagonal eigen- 
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[A]{y> = 


11 


22 


33 


m 


m 


mm 


{y} = A{y} (27) 


and the matrix coefficients are theoretically given by the simplified recurrence 
formulas 


a i,i = {v i }T[ B^v i } 

{ W = [B](v.} - a.^} - d.^} U-1,1 


d i+l = [ ^i+i > T 1 { v ±+i > 3 1 7 2 


(28a) 


* V i+l* d... * V i+l^ ; 
1+1 


1=1, m-1 


(28b) 


where the sequence is initialized by choosing a random starting vector for {v^} 
and setting d^ = 0, {v q } = {0} • 

Although the final off-diagonal term, given by Equations (28a) is 

not needed to construct the tri-diagonal reduced matrix, it is calculated and 
saved for use in establishing error bounds (see Section 6). In addition, the 
above algorithm is modified in the computer program as follows: 

1. Each vector {v.,-.}, calculated by Equations (28b) is reorthogonalized 
to all previously computed {v}-vectors and eigenvectors, as des- 
cribed in Section 4, before re-entering Equations (28a). 
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2 


. The actual size, m, of the reduced eigenproblem is established by 
the criteria given in Section 2.3, one of the restrictions being 
that it cannot exceed the rank of the matrix [M aa ] for vibration 
mode problems or [K^ a ] for buckling problems. 


The eigenvalues. A, and eigenvectors, {y} of Equation (27) are extracted 
using a Q-R algorithm and eigenvector computational procedure similar to that 
described in the NASTRAN Theoretical Manual (Reference [13]). They are then con- 
verted to physical form as follows: 


X 


i 



(buckling problems) (29a) 


-2 


a). 


i 



(unshifted vibration (29b) 

mode problems) 


-2 

to. 

X 



(shifted vibration 
mode problems) 


(29c) 


{<}>.}= [C "*"] T [V]{y.} (buckling or unshifted (29d) 

vibration mode problems) 

= [V]{y^} (shifted vibration (29e) 

mode problems) 


T 

In addition, the vectors {y^} are normalized ({y^} {y^} = 1) f° r convenience in 
establishing error bounds (see Section 6) . 

2.3 Criteria for the Size of the Reduced Eigenvalue Problem 

The total number of eigensolutions, including any existing rigid body 

modes, is equal to the rank, r, of the matrix [M] = [M ] or [M] = [K^ ], depend- 

aa aa 

ing on whether a buckling or vibration mode problem is being solved. Thus, the 
size of the reduced problem, m, cannot be greater than r. If, for example, [M] 
is diagonal, then the maximum permissible value of m is equal to the total number 
of non-zero diagonal entries. In addition, if f eigensolutions have previously 
been computed by NASTRAN, these must be swept out of the problem by making all 
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{v} vectors orthogonal to the previously computed eigenvectors. This implies 
that the maximum size of the reduced problem is further reduced to 

r = r - f (30) 

As a result of numerical experiments and application experiences 
(References [8, 9, 14]), it has been found that in cases where m « r, a first 
grouping of more than ra/2 eigenvalues closest to the shift point are in 
accurate agreement with the corresponding number of exact eigenvalues. The 
remaining reduced-system roots are spread across the remaining exact eigen- 
spectrum. 

Thus, if the user requests a total of q eigenvalues closest to a speci- 
fied numerical value*, the order of the reduced problem solved by NASTRAN will be 

m = min[2q +10, r] (31a) 

where 

q = q - f (31b) 

It should be noted that in all cases m £ r, and whenever m is set 
equal to r, all the eigensolutions of the unreduced problem, which have not been 
previously computed, are generated. 

From the above discussion it is apparent that, depending on the number of 
eigensolutions requested, the rank of [M] may have an influence on the order of the 
reduced problem which is generated by the program. Although the user should never 
request more than r eigensolutions, the value of r may not always be a simple 
matter to calculate, particularly in buckling problems. Therefore, the computa- 
tional scheme has been designed to resolve the question of rank in the following 
manner: 


*q includes the number of previously computed eigensolutions, f. These consist 
of modes generated prior to a restart plus rigid body modes generated by using 
a SUP0RT card in the bulk data deck. 
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1. The matrix [M] is first checked for inordinately small off-diagonal 
elements, i.e., those which lie in the round-off range of the com- 
puter. These terms can sometimes introduce artifically large physi- 
cal eigenvalues and are therefore eliminated. Any off-diagonal 
element for which 


-2t 



is set equal to zero, where t is the number of decimal digits 
carried in the computation of the [M] matrix. 

2. The number of non-null columns or rows of the above modified [M] 
matrix is counted and designated as n. Since the rank, r, cannot 
be greater than n, the program initially sets r = n - f in Equa- 
tions (31) to provide a tentative size, m, of the reduced eigenvalue 
problem. If the user has asked for more than r eigensolutions , he 
is given a message that too many eigensolutions have been requested 
and that the program will try to find all the existing solutions. 

(See Section 9 for a complete list of user messages and diagnostics) . 

3. If the reorthogonalization tests (see Section 4) fail for some 

vector this is an indication that a null vector has been 

generated because the maximum number of r-f linearly independent 
{v} vectors have already been obtained. The recurrence sequence 
is then terminated and the order of the reduced eigenproblem is 
further reduced to m = i. 
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3. 


CHOICE OF THE INITIAL TRIAL VECTOR AND RESTART VECTORS 


Prior to tridiagonal reduction, the original eigenvalue problem 


is cast in the inverse form 


or 


where 


[K - A M]{(j)} = 0 

3 . 


[B] {X} = A[D] {X} 
[B] {X} = A{X} 



( 1 ) 

( 2 ) 

(3) 

( 4 ) 


From Equation (1) it can be seen that if the problem size is designated 
as n and r is the rank of the [M] matrix, then there are (n-r) spurious 
eigensolutions corresponding to A = 0 or, equivalently, A^ -* . If, for example, 

a relatively small number of nodal masses or differential stiffness elements are 
employed, then [B] contains a large multiplicity of zero eigenvalues which are of 
no interest and cause numerical difficulties. These are eliminated from the re- 
duced tridiagonal problem by employing a constrained sub-set of v-vectors having 
zero projection on the set of eigenvectors associated with A = 0. 

To accomplish this, use is made of the fact that any non-null vector {v^}, 
generated from any other non-null vector {w}, through 


{v^} = [B] {w} (5) 

will contain no components of the eigenvectors corresponding to A = 0. This can 
be seen as follows. 

First, express {w} as a linear combination of all the eigenvectors of [B] : 


- 17 - 


( 6 ) 


{w} = I c {X } + l c {X } 
i=l 1 i=r+l 


where {X..}, (i = 1, r) , are the eigenvectors for A ^ 0 and {X^}, (i = r + 1, n) , 
are the eigenvectors for A = 0. 

Next, substitute Equation (6) into Equation (5), giving 


{v,} - I c . [ b] {x. } + l C [B]{x :} 
i=l i=r+l 


(7) 


Since 


and 


[B] {X. } = A. {X. } ; i = 1, r 

(8) 

1 XI 

[B] {x^} =0 ; i = r+1, n 

(9) 


it follows that 


'V ■ £ «iW 


( 10 ) 


Thus, {v^} contains no eigenvectors corresponding to A - 0, and is a null vector 
only if {w} happens to be a A = 0 eigenvector. Further, it is easy to see that 


after {v^} is normalized 


{v } = [ }, (11) 

1 {v 1 } T [D](v 1 > 1 

the next vector generated, 

{v 2 > = [B ] > - a 1 2^}, (12) 

as well as all subsequent trial vectors will be free of A = 0 eigenvectors. 

Employing a somewhat similar argument, it can be shown that the most de- 
sirable initial trial vector, {v^}, and hence the vector {w} from which it is 
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generated, should contain all components of the eigenvectors for A ^ 0. How- 
ever, since there is generally no a priori knowledge of the modal matrix, {w} 
should be selected in such a way so as to make it as "irregular" as possible with 
respect to the system of eigenvectors so that it is most likely to contain a mix- 
ture of all the mode shapes. It has been found that this is best achieved by 
using a random or pseudo-random number generator to obtain the elements of 
{w}. 


If the vector {w} selected in the above manner should, by some chance, 

be deficient in eigenvector components, then a null v-vector may be generated 

at some point by the recurrence algorithm (Equations (28), Section 2.2). In the 

context of a finite digit computer, this is indicated by the appearance of an 

off-diagonal term, d i+ ^, which is exceedingly small compared to the corresponding 

diagonal term, a. in the reduced tridiagonal matrix. The test used is that 
x , i 

U i+1 l iio 2 ' 1 |a 1>± l (13) 

implies that the newly generated vector {v i+ ^}, is null, where t is the number 
of decimal digits carried by the computer. In this event, d^ + ^ is set exactly 
equal to zero and a new restart vector is employed for {v^ + ^}. This vector is 
generated exactly as in the case of the initial trial vector {v^}, but using a 
different pseudo-random number seed. The recurrence algorithm for generation of 
the v-vectors is subsequently continued in the usual manner until the required 
number of vectors has been generated. 

A further constraint on each v-vector, including the initial one, is that 
it be orthogonal to all previously generated v-vectors and previously calculated 
eigenvectors*. The imposition of these additional constraints is discussed in 
the next section. 


*Previously calculated eigenvectors may be available from the following sources: 

(1) The specification of fictitious free body supports on a SUP0RT card, which 
causes an equal number of rigid body modes to be automatically generated by 
NASTRAN prior to entering an eigenvalue extraction routine. 

(2) Checkpointing of previously obtained eigensolutions followed by a restart 
to obtain additional eigensolutions. 
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4. SW EEPING- OUT OF PRE VIOUSLY OBTAINED EIGENVECTORS AND REORTHOGONALIZA- 

TiON OF THE TRIAL VECTORS 

Assume that a combination of f rigid body and non-rigid body eigen- 
vectors have already been extracted prior to the current application of the 
tridiagonal reduction method. Let these vectors be designated by {xp, {xp. 

In order to avoid regenerating these previous eigensolutions, which 
would be inefficient, the initial trial vector, {v^}, obtained in the manner just 
described, should be made orthogonal or "swept clean" of these eigenvector com- 
ponents. This is theoretically accomplished by setting 


{vp = {v 


,} - l l {X. } T [D j {v. } ] {X. } 
j=l 3 X 3 


( 1 ) 


where is the swept initial vector to be used in place of {v^} following its 

normalization. It can then be shown that all succeeding v-vectors generated by 
the recurrence algorithm (Equations 28, Section 2.2) form a theoretically ortho- 
gonal set which does not contain components of {X.}; j = 1, f. However, it has 
been shown (Reference [15]) that the v-vectors degrade rapidly as the co input at ions 
proceed, such that the later vectors are far removed from orthogonality to the 
earlier ones. This is caused by unavoidable computational round-off which, because 
of repeated multiplications by the unreduced eigenmatrix, [B], tends to amplify 
the contributions of the eigenvector components nearest the shift point in the 
calculated trial vectors. Thus, unless sufficiently accurate orthogonality of the 
trial vectors is maintained, they will be excessively rich in the modes near the 
shift point, and the solution of the reduced tridiagonal eigenvalue problem will 
yield a false bunching of eigenvalues around this point in the eigenspectrum. 


To correct this problem, Gregory (Reference [6]) experimented with the use 
of higher-precision computer operations, but found only marginal improvement in 
the final results. Later, Lanczos suggested a reorthogonalization of the type 


{v i+ i } - {v ±+i ) - H V 


[D]{v )Hv > - l 

3 j=l 


[{XjKfDHv^mx.} (2) 
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where {v i+1 } is calculated by the unmodified recurrence algorithm and {v i+ ^} is 
an improved vector- While this improves matters substantially, it still does not 
eliminate the precision problem adequately. However, Ojalvo and Newman (Reference 
[8]) found that the introduction of an iterative reorthogonalization loop can make 
the trial vectors as orthogonal as necessary for extremely large systems. The 
procedure is as follows: 


The vector {v.,..}, obtained from either the recurrence algorithm or the 

fO} 

pseudo -random number generator (see Section 3) is denoted as {v^ + £} and reortho- 
gonalized with respect to all the previously obtained vectors. This is accom- 
plished by iterating. 


<v£> 


w i+l J 



[{X_.} T [Di{v^}]{X.} 

[{X.} T [D]{v^}]{X.} 


(3) 


until an acceptable vector 


{v 


(s+1) 

i+1 


} = 




[ IX. } T [D] (v^ } ] {Xj > 


(4) 


is found which satisfies the orthogonality criteria 

taax | {v. } T [D] {v^ } | £ 10 2 (5a) 

1 < j < i 3 

max | {X. } T [D] {v£f} } | < 10 2_t (5b) 

1 < j < f 3 

where t is the number of decimal digits carried by the computer. If, for some 
vector, the above criteria are not satisfied after a set number of iterations , then the 
program assumes that a new trial vector cannot be generated, and a reduced eigen- 
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value problem of order m = i is solved, as discussed in Section 2.3. 

If the above criteria are met, then the resulting vector is normalized 
and set equal to the new normalized trial vector, i.e.. 


( 6 ) 


This new vector is used to compute the next off-diagonal term in the reduced tri- 
diagonal matrix from the formula 

d i+l = <v. +1 } T [BHv i } . (7) 

However , if 

K +1 l iio 2_t K.ii. «> 

it is probable that {v^ + ^} iS a nul -*- vector, possibly because the maximum number 
of linearly independent vectors, corresponding to the rank of the problem, has 
been exceeded. In this event a reduced eigenvalue problem of order m = i is 
solved , as above . 

If the criteria given by Equations (5) and (8) are all met, then the new 
normalized vector {v^_^} is used to continue the reduction algorithm. 
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I 


5. CRITERION FOR THE NEGATIVE SHIFT PARAMETER, g 2 

Before tridiagonal reduction, the eigenvalue problem for natural fre- 
quencies in the neighborhood of zero is 


where 


and 


[b] {x} = A{x} 

(1) 

[B] = [C _1 ][M 1[C~ 1 ] T 

da 

(2) 

{X} = [C] T {<J>} 

(3) 

IK «. + “\a 1 - tcltcl 1 

(4) 

A ■ ..2 t 2 

(5) 


The criterion involved in the choice of the shift parameter, a , is that 
it be large enough to render a possibly singular stiffness matrix non— singular (to 
the extent that Cholesky symmetric decomposiion can be performed accurately in a 
finite digit computer) , and small enough to prevent troublesome clustering of the 
eigenvalues, A^. As an approach toward solving this problem, it is helpful to 
note that, given a symmetric positive-definite matrix [J], we really obtain 


[C][C] T = [J] + [6J] 


( 6 ) 


when it is factored, because of computer rounding errors. According to Wilkinson 
(Reference [16]), the following inequality is almost always satisfied, 

llfijlli < n(10 1_t ) I I j I I x (7) 

where the above Holder-one matrix norms are equal to the maximum row sums of 
absolute values, n is the order of the matrix, and t is the number of decimal 
digits carried by the computer. 
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Thus, in order to render the matrix [K + a M ] non-singular, even 

aa aa 2 

when [K ] itself is singular, some matrix norm of the modification a [M ] 
aa aa 

should be appreciably larger than the corresponding norm of [6 k ], where the 

aa 

latter matrix is in the "noise level" of the computer. On the basis of Equa- 

2 

tion (7), the minimum a should therefore satisfy the requirement 


af . | |M | | > n(10 1 fc ) | |K || 
mm 1 ' aa 1 1 1 — x 7 11 aa M 1 


or 


a 2 . > 

mxn — 


f" max l |K 1 1 

n(10 l -‘) J' - 1 ij 

|~ “ax l |m I] 

Ll 1 i < n j=l 1J J 


where K. . and m. . are the elements of [K ] and [M ]. 
1 3 aa J aa 


( 8 ) 


(9) 


To further enhance the removal of possible singularities in the stiffness 

matrix, a factor of ten is applied to the right side and the maximum ratio of 

? 

diagonal matrix elements is used. The resulting value of oT . is 


a 2 . = n(10 2 t ) |-ii| ; 

mxn ‘m . 'max 


xx 


m ±i * 0 


( 10 ) 


A final requirement imposed on the shift parameter is that it be large 
enough to effect alterations in the last two-thirds of the significant digits in 
each diagonal term of the unmodified stiffness matrix or, equivalently. 


a 2 = 10 |— | . ; 

o 'm. . 'mm 

li 


m. . ^ 0 

li 


(ID 


Consequently, the actual value of the shift used by the program is 

2 2 2 
a = max (a . , a ) 
mm o 


( 12 ) 
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If the resulting modified stiffness matrix is still singular, as Indicated 

2 

by a failure of the Cholesky decomposition process, the above value of a is 
multiplied by a factor of 100, a maximum of two times in an attempt to render 
the undecomposed matrix non-singular. If this procedure fails, the problem exe- 
cution is aborted and the user is informed that the singularities cannot be re- 
moved from the stiffness matrix. 
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6 . 


ERROR BOUNDS ON THE COMPUTED EIGENVALUES 


Once the inodes of the reduced problem have been extracted, close upper 
and lower bounds on the eigenvalue errors can be obtained rather economically. 
This provides the user with an a posteriori check on the number of accurately 
calculated eigenvalues which, in most cases, will be greater than the number re- 
quested. 

In carrying out the development, it is convenient to write the eigen- 
problem in its inverse form 


where 


[B] {x} = A{X} 


( 1 ) 



( 2 ) 


Considering some approximation, {X^}, to an exact eigenvector, {X^}, it 
can be expanded in terras of the exact modal matrix as follows: 


where 


and 


{x ± } = [Z3 [b ± l 


( 3 ) 


[z] = t tx 1 > ; 

nxn 


{x 2 > i 


: {x n >] 


( 4 ) 



( 5 ) 


is a vector of scalar coefficients. 

If the eigenvalue approximation associated with {X_^} is denoted by A^, 
then the residual vector for the ith modal approximation is given by 
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{R ± > = [B - A i [I]]{X i > 

= tB - A 1 [I]][Z]{b 1 } 


However 

[B][Z] = [Z]hA d J 


where ["A^] is the diagonal matrix of eigenvalues, A^, i = 1, n. 
It follows that 


( 6 ) 


( 7 ) 


= [Z][['A d J - A i [I]]{b.} (8) 

Denoting a weighted root-mean-square residual by 

r i = ({ R i > T [ D ]tR i » 1/2 , (9) 

it can be seen that 

= {b.} T [rA d J - A 1 [l3][Z] T [DJ[Z][rA d ^] - A ± [I]]{ b ± } (10) 

However, without any loss of generality, the eigenvectors can be orthonormalized 
so that 

{x j } T [D]{x i > = { (id 


The modal matrix is therefore unitary, i.e.. 


and 


[Z] T [D][Z] = [I] 


2 
r . 
i 


= {b.} T [[~A ] - A [I]] 2 {b } 

1 a 1 1 


( 12 ) 


n 


= l 

j=l 


b 2 .. (A. - A.) 2 
J i 


( 13 ) 
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where {b J4 } refers to the ith element of the vector {b.}. 
ij 1 

At this point it should be noted that NASTRAN also normalizes the eigen- 
vectors {y.} of the reduced tridiagonal problem, so that 


fw iT r 


where 


Thus, 


and it follows that 


{X.} T [ 


or 


X . - X 
2.) 

ai o 


[v] t [d] [ v]{y ± } = l 

(14) 

,} i..J {v >] 

m 

(15) 

n 9 

i b?. . i 
-i ^ 

(16) 

(A. - ip 2 

(17) 

eigenvalue , A^, it can be 

seen that 

£ |r.| 

(18) 

1| 1 l-l 

(19) 


A. 

x 


Since X . is close to X . when A. is close to A., Equation (19) 
ax ai x x* 


implies that 


X . r. 

|i - x^l < b — 

ai A.(l + X A.) 

X ox 


( 20 ) 


which is a measure of the maximum relative error in the ith physical eigenvalue*. 


The residual, r^, in the above equation can be evaluated quite easily via 


the following approach. 


*This test is obviously invalid when X a ^ ~ 0, i.e., a rigid body solution. In 
this event, the computed value of X ai is, in itself, a measure of the absolute 
error in the physical eigenvalue and no further accuracy information is needed. 
The existence of an eigenvalue as an additional rigid body mode not requested 
on a SUP0RT card in the NASTRAN Bulk Data Deck is detected by the criterion 
|X ai | < , in which case Equation (20) is bypassed. 


- 31 - 



Set 


{X.} = [Vl{y ± > 


in Equation (6) . Then 


{R.} = [ [B ] [V] - A ± [V]]{y i }. 

However, from Equation (26) of Section 2.2, it is seen that 


[B][V] - [V] [AJ + 


where [A] is the tridiagonal reduced matrix and 


I<W - 1{0 > ! {0} <Vi )] 


Therefore, 


But 


so that 


and 


(V - [ [v] IA] + a^IG^l - Ajvnty,} 


[A]^} = A i {y i J, 


{R i> ■ WVi Hy i } - 


r i • < K l> T t D l« R i} - 4n 


It can be seen that 




0 

0 


0 


0 0 

0 0 

• • 

0 1 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 
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so that finally 


where y . 

mi 

Therefore, 


2 2 2 
r i = Vl y mi » 


is the last element of the vector {y_^}. 
Equation (20) assumes the form 


( 30 ) 


A . I d * y - 1 

ai | < 1 nrfl mi 1 

X ai “ 1 A. (1 + A A.) I 
1 1 ox' 


(31) 


Thus, it is seen that the eigenvalue errors are proportional to d^^, which is 
the next off-diagonal term that would be generated, had the reduced tridiagonal 
matrix [A] been increased from order m to order m+1. Equation (31) shows that 
this term is further modified by a weighting factor y which is the last term 
in the reduced-system eigenvector associated with A^. 

The use of the above error bound formula as a criterion for selecting 
acceptable eigensolutions is described in Section 7. 
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7. 


SUMMARY OF COMPUTATIONAL PROCEDURES AND PLOW CHARTS 


Flow diagrams illustrating the computational procedures are shown in 
Figures 1 and 2. The details of each block are summarized below. 

2 

(1) Calculate Small Negative Shift Parameters, a , (See Section 5) 

In the case of unshifted vibration mode problems the negative shift para- 
meter for removing possible singularities is found from: 


where 


and 



( 1 ) 

(2a) 


(2b) 


K. . and m.. are the diagonal elements of [K ] and [M ], respectively, n is 
n li aa aa 

the number of {u } degrees of freedom, and t is the number of decimal digits 
a 

carried by the computer. 

(2) Zero-Out Excessively Small Elements of [M] Matrix (See Section 2.3) 

a. Compare the magnitudes of all off-diagonal elements of [M] with 
the corresponding diagonal elements to determine whether 


b^-l < 10 2t/3 ; i = J, m 4 0 . (3) 

m. . — li 

n 

b. Set m _ =0.0 for every off-diagonal element satisfying the above 
criterion. 

(3) Establish Tentative Reduced Problem Size (See Section 2.3) 

a. Count the number, n, of non-null columns or rows in the above modi- 
fied [M] matrix and set 
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(4) 


r " n - f 

where f is the number of previously computed eigensolutions . 

(b) Calculate a tentative size, m, of the reduced eigenproblem from 


m = min[2q + 10, r] (5 a ) 

where 

q = q - f (5b) 

and q is the total number of accurate eigenvalues requested by the user, 
including previously computed modes. If q > r, warn user that too many 
eigensolutions have been requested - program will try to find all existing 
solutions. 

(4) Construct Factors of [K] Matrix (See Section 2.1) 

a. Set 



(i) 

[K] = [K - W* M l 
aa o aa 

(Shifted Vibration 
Mode Problems) 

(6a) 

or 

(ii) 

KI - [K aa + Aj 

(Unshifted Vibration 
Mode Problems) 

(6b) 

or 

(iii) 

[K] - [K aa ] 

(Buckling Problems) 

(6c) 

b. 

Perform a 

non-square root decomposition: 




[K] = [L] [" d 

T 

-HL] 1 

(7a) 

for 

case (i). 

or a Cholesky symmetric 

decomposition : 
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(7b) 


[K] = [C][C] T 

for cases (ii) and (iii), using real arithmetic without pivoting. Save 

the triangular and diagonal factors. If the decomposition for case (i) 

fails or the decompositions for cases (ii) and (iii) fail after two in- 
2 

creases in a by factors of one hundred, then the program is aborted 
and a fatal error message is issued, indicating that there are unremov- 
able stiffness matrix singularities. 

(5) Execute Tridiagdnal Reduction Algorithm (See flow diagram for this block. 

Figure 2) 

(5.1) Initialize the Recurrence Algorithm (See Section 2.2) 

Initialize the vector index to i = 0 and set 


{v Q } = {0} (8) 

where {v } is an (nxl) null vector, 
o 

(5.2) Generate a Starting or Restart Vector and set d i+ ^ = ^.0 
(See Section 3) 

(a) Construct an n-element vector {w} using a pseudo-random 
number generator. 

(b) Solve for an un-normalized trial vector from the equation 


where 


or 


or 


{v ±+i } = tB]{w} (9) 

[B] = [L" 1 ] T hd.]" 1 [L" 1 ]tM aa ] (easel) (10a) 

a.a 

[B] = [C 1 ] [M HC-V (case ii) (10b) 

aa 

[B] = [C -1 ][K d ] [C (case iii) (10c) 
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Forward and backward passes are used to perform the above inverse 
operations . 


(c) Normalize the above vector: 


MS 1 


{v l+l } [D]{v 1+1 > 


1 1/2 <W 


( 11 ) 


where 


ID] = [M l 
aa 

(case i) 

(12a) 

[D] = [I] 

(cases ii and iii) 

(12b) 


(d) Set " 0*0 and Proceed to block (5*5). 

(5.3) Create One Approximate Trial Vector and One Diagonal Coefficient 
(See Section 2.2) 

The recurrence algorithm is: 


3 i,i = 


(13) 


{v i+ i> = - d.tv.^} 


(14) 


\+l = Hv i+ i> T [ D Hv i+1 }] 1/2 


(15) 


trill > = — <v_. „ > 


i+1 


i+1 


* 1+1 


where 


[B] = [D][B] 


( 16 ) 


(17) 


and {vf^} I s an approximation to the new trial vector. 
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(5.4) 


First Normalization Test (See Section 3) 
The test is 


• d i+il > 1q2 l a i,il < 18 > 

Pass: Proceed directly to block (5.5) 

Fail: Return to block (5.2), generate a new restart vector for } , 

and then proceed to block (5.5). 


(5.5) Iterate to Obtain Orthogonalized Vector (See Section 4) 


Designate {X_.} , j 
eigenvectors. Perform the 


= 1, f as previously calculated and stored 
iterations. 


{v m 1)> ■ (v ^ ) - 1 


1+1 A J 


i+1 j 


f 

- I 


j=1 IX j } T [D]{v^}]{X j }; s = 0,1,2. 


until 


and 


or 


max |{v > T [D]{v^}| £ 10 2 C 

1 < j < i 3 


max |{X } X [D]{v^}| < 10 2_t 

1 < j < f 


(19) 


( 20 ) 


s = 14. 


If the orthogonality criterion. Equation (20), is satisfied, proceed to 
block (5.6). Otherwise, set the problem size, m, equal to i and pro- 
ceed to Exit. 

(5.6) Normalize the Orthogonalized Trial Vector (See Section 4) 
Compute 
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( 21 ) 


(v 

i+1 


} = 


{ v ( s+1 )} 

w i+l 


[<v lS 1)}lD){v -^ 1)}ll/2 


i+1 


This is the new orthogonalized and normalized trial vector. 

(5.7) Second Normalization Test and Creation of Off-Diagonal Coefficient 
(See Section 4) 

(a) Compute the next off-diagonal term in the reduced tridiagonal 
matrix from 


d i+i = {v ±+i } [B]{v ± } 


(b) Verify whether the following test is met: 


> i ° 2 ' t |a i,i 


( 22 ) 


(23) 


If it has, set i = i+1 and return to block (5.3) for continuation of 
the recurrence algorithm. If the test fails, set m = i to reduce the 
problem size and proceed to Exit. Issue message to user that only i 
modes can be obtained - more than r - f modes may have been requested. 

(6) Solve Reduced-System Eigenvalue Problem (See Section 2.2) 

a. The coefficients a_ - , a 00 , ... , a and d 0 , d^, ...., d , 

11 22 mm 2 3 m 

computed in block (5) are interpreted as the following symmetric, tridiagonal 
array: 


a ll d 2 


22 


[A] = 


33 


(24) 


m 


tn 


mm 
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b. The mth order eigenvalue problem 

[Aj{y} = A{y} (25) 

is solved for the eigenvalues, and eigenvectors {y^} using a Q-R algorithm 

and eigenvector computational procedure similar to that described in the NASTRAN 
Theoretical Manual (Reference [13]). 

c. The reduced system eigenvectors are normalized so that 

{y i }T{y i } = 1; 1 = m (26) 

(7) Compute Maximum Eigenvalue Errors (See Section 6) 

a. The maximum absolute relative errors in the computed physical 

eigenvalues are obtained from 


?i 


I (W (y mi ) 1 . 
l A i (1 + A o A ± ) | 


i = 1, m 


(27a) 


where d , - is the last off-diagonal term computed in block (5,3) and y . is 
tn+l w mi 

the last element in the vector {y^}. If the physical eigenvalue 1/A^ + X q> 
corresponds to a rigid body mode, the above computation is invalid and therefore 
bypassed. A rigid body mode is assumed to occur whenever 

— *- ! ^. 10_t/3 (27b) 

and is denoted by setting the relative error, equal to a flat zero. 

b. The eigenvalues are processed in order of increasing distance 

from the center of range of interest, X q , to determine whether their associated 
values meet an acceptable relative error tolerance set by the user on the 
EIGR or EIGB bulk data card (the default value is .001/n where n is the order 
of the stiffness matrix). The first eignevalue not meeting the tolerance test. 
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as well as all subsequent eigenvalues further removed from the center of interest, 
are considered to lack sufficient accuracy and are therefore rejected. 

c. Acceptance eigenvalues obtained in the above manner are reordered 

in terms of increasing physical value for subsequent processing by NASTRAN. 

8. Compute Physical Eigenvalues and Eigenvectors (See Section 2.2) 


The mathematical eigenvalues, 
to physical form as follows: 



(4^} = [C _1 ] T [V]{y i } 
{* ± } = [V]{y i > 

where 

[V] = [W x 


and eigenvectors, {y.^}, are converted 


(Buckling Problems) (28a) 

(Unshifted Vibration (28b) 

Mode Problems) 

(Shifted Vibration (28c) 

Mode Problems) 

(Buckling or Unshifted (29a) 

Vibration Mode Problems) 

(Shifted Vibration (29b) 

Mode Problems) 

2 > {%}]* (30) 
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I 




Figure 1. Overall flow diagram for tridiagonal reduction method 
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5.3 Create one 

approximate trial 
vector and one 
diagonal coeff. 


Initialize the 
Recurrence 
Algorithm; i = 0 


^ i = 
Problem 


5.2 

Generate a starting 

or restart vector 

and set d. ... = 0.0 
l+l 

l 

5.5 

Iterate to obtain 

orthogonalized 

vector 


First \ 
.normalization, 
X test 


Unsatisfied 


Orthogonality^ 
^criterion / 


Reduce 
problem 
size to i 


yp Satisfied 

Normalize the 
orthogonal ized 
trial vector 


Second 

normalization 
test and creation 
of off-diagonal 
coeff. 


i = i+1. 


Figure 2. Flow diagram for block 5, execute tridiagonal reduction algorithm 






8 . 


NASTRAN USER’S INSTRUCTIONS 


The following pages show modifications of the EIGR and EIGB cards in 
the NASTRAN Bulk Data deck which accommodate user implementation of the tridia- 
gonal reduction method for real eigenvalue analysis. The modifications are 
constituted of additions to the standard user instructions and are underscored 
for ease in identification. 

When the tridiagonal reduction method is invoked, the F2 or L2 parameter 
on these cards represents the maximum allowable value of the computed relative 
error in a physical eigenvalue (see Section 6). If this value is exceeded, the 
associated eigensolution is not accepted for further processing by NASTRAN. A 
detailed list of the maximum relative errors in the computed eigenvalues can be 
obtained by requesting DIAG 16 in the NASTRAN Executive Control Deck (see Section 
9). 
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BULK DATA DECK 


Input Data Card EIGB Buckling Analysis Data 

Description : Defines data needed to perform buckling analysis. 

Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

EIGB 

SID 

METHOD 

LI 

L2 

NEP 

NDP 

NDN 

E 

+abc 

EIGB 

13 

DET 

0.1 

2.5 

2 

i 

1 

0.0 

ABC 


+abc 

N0RM 

G 

C 







+BC 

MAX 










Field 


Contents 


SID 

METH0D 


LI, L2 


NEP 


NDP, NDN 


Set identification number (Unique integer > 0) 

Method of eigenvalue extraction, one of the BCD values "INV", 
"DET", "FEER" , "UINV", or "UDET" 

INV - Inverse power method, symmetric matrix operations 

DET - Determinant method, symmetric matrix operations 

FEER - Tridiagonal reduction method, symmetric matrix 
operations 

UINV - Inverse power method, unsymmetric matrix operations 

UDET - Determinant method, unsymmetric matrix operations 

Eigenvalue range of interest (Real; LI < L2 > 0.0) For METH0D 
= "FEER", LI is ignored and LZ is the acceptable relative error 
tolerance on eigenvalues, in percent (Default is .001/n where n is 
the order of the stiffness matrix) (Real > 0.0) 

Estimate of number of roots in positive range. Desired number 
of eigenvalues of smallest magnitude for METH0D = "PEER 11 
(Default is automatically calculated to extract at lease one 
accurate mode) (Integer > 0) 

Desired number of positive and negative roots (Default = 3 NEP) 
(Integer >0). Ignored for METH0D = "FEER" 
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E 


Convergence criteria (optional) (Real > 0.0) 


N0RM 

G 

C 

Remarks 

1 . 

2 . 

3. 

4. 

5 . 

6 . 

7. 


Method for normalizing eigenvectors, one of the BCD values "MAX" 
or "P0INT" 

MAX - Normalize to unit value of the largest component in the 
analysis set 

P0INT - Normalize to unit value of the component defined in 

fields 3 and 4 defaults to "MAX" if defined component, 
is zero 

Grid or scalar point identification number (Integer > 0) (Re- 
quired if and only if N0RM « "P0INT") 

Component number (One of the integers 1-6) (Required if and 
only if N0RM = "P0INT" and G is a geometric grid point) 


Buckling analysis root extraction data sets must be selected in the 
Case Control Deck (METH0D = SID) to be used by NASTRAN. 


The quantities LI and L2 are dimensionless and specify a range in 
which the eigenvalues are to be found. An eigenvalue is a factor by 
which the prebuckling state of stress (first subcase) is multiplied to 
produce buckling . If METH0D = "FEER", LI is ignored and L2 repre- 
sents the maximum upper bound, in percent, on A gxACT ~ ^ or 

acceptance of a computed eigensolution . 


The continuation card is required. 


See Sections 10.3.6 and 10.4.2.2 of the Theoretical Manual for a dis- 
cussion of convergence criteria. 

If METH0D = "DET", LI must be greater than or equal to 0.0. 

If N0RM = "MAX", components that are not in the analysis set may have 
values larger than unity. 

If N0RM = "P0INT", the selected component must be in the analysis set. 
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BULK DATA DECK 


Input Data Card EIGR 

Description: Defines data needed to perform real eigenvalue analysis. 

Format and Example: 


EIGR 

SID 

METH0D 

FI 

F2 

[ NE 

ND 

NZ 

E 

1 -fabc 

EIGR 

13 

DET 

1.9 

15.6 

10 

12 

o 

r i.-3 

! ABC 


+abc 

N0RM 

G 

c i 







+BC 

P0INT 

32 

4 








Field Content 

SID Set identification number (Unique integer > 0) 

METH0D Method of eigenvalue extraction, one of the BCD values "INV", 

"DET", "GIV", "FEER" , "UINV", or "UDET". 

INV Inverse power method, symmetric matrix operations. 

DET Determinant method, symmetric matrix operations, 

GIV Givens method of tridiagonalization. 

FEER Tridiagonal reduction method, symmetric matrix 
operations. 

UINV Inverse power method, unsymmetric matrix operations. 

UDET Determinant method, unsymmetric matrix operations. 

FI, F2 Frequency range of interest (Required for METH0D = "DET", 

"INV", "UDET", or "UINV") (Real > 0.0; FI < F2) . Frequency 
range over which eigenvectors are desired for METH0D = "GIV". 
The frequency range is ignored if ND > 0, in which case the 
eigenvectors for the first ND positive roots are found. (Real 
FI < F2). If METH0D = "FEER", FI is the center of range of 
interest (Default is FI = 0.0) (Real _> 0.0), and F2 is the 
acceptable relative error tolerance on f requency-squared , in 
percent (Default is .001/n where n is th^ order ojf the stiffness 
matrix) (Real > 0.0) 

NE 

Estimate of number of roots in range (Required for METH0D = 
"DET", "INV", "UDET", or "UINV", ignored for METH0D = "FEER") 
(Integer > 0) 
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ND Desired number of roots for METH0D = "DET", "I NV", "UDET" , 

or "UINV", (Default is 3 NE) (Integer > 0). Desired number 
of eigenvectors for METH0D • - "GIV" (Default is zero) (Integer 
0). Desired number of roots and eigenvectors for METH0D = 
"FEER" (Default, is automatically calculated to extract at least 
one accurate mode) (Integer > 0) 

NZ Number of free body modes (Optional - used only if METH0D = 

"DET 11 or "UDET") (Integer > 0) 

E Mass orthogonality test parameter (Default is 0*0 which means 

no test will be made) (Real _> 0.0). 

N0RM Method for normalizing eigenvectors, one of the BCD values 

"MASS", "MAX" or "P0INT" 

MASS - Normalize to unit value of the generalized mass 

MAX - Normalize to unit value of the largest component in 

the analysis set 

P0INT - Normalize to unit value of the component defined in 

fields 3 and 4 - defaults to "MAX" if defined com- 

ponent is zero 

G Grid or scalar point identification number (Required if and only 

if N0RM = "P0INT") (Integer _> 0) 

C Component number (One of the integers 1-6) (Required if and 

only if N0RM «= "P0INT" and G is a geometric grid point) 


Remarks : 


1 * 

2 . 

3. 

4 * 


Real eigenvalue extraction data sets must be selected in the Case 
Control Deck (METH0D = SID) to be used by NASTRAN. 

The units of FI and F2 are cycles per unit time. If METH0D = 
"FEER", F2 represents the maximum upper bound, in percent, on 

- - j- *'■'» ----- — i — - — r : — : 

I^FEER— 


- 1| for acceptance of a computed eigensolution . 


The continuation card is required. 

If METH0D = "GIV", all eigenvalues are found. 


5. If METH0D = "GIV", the mass matrix for the analysis set must be 

positive definite. This means that all degrees of freedom, including 
rotations, must have mass properties. 0MIT cards may be used to re- 
move massless degrees of freedom. 


- 49 - 


6 . 


A nonzero value of E in field 9 also modifies the convergence 
criteria* See Sections 10.3.6 and 10.4*2.2 of the Theoretical Manual 
for a discussion of convergence criteria. 

7. If N0RM “ "MAX", components that are not in the analysis set may have 
values larger than unity. 

8. If N0RM = "P0INT", the selected component must be in the analysis set. 

9. If METH0D = "GIV" and rigid body modes are present, FI should be 
set to a small negative number rather than zero if the rigid body 
eigenvectors are desired. 

10. The desired number of roots (ND) includes all roots previously found, 
such as rigid body modes determined with the use of the SUP0RT card, 
or the number of roots found on the previous run when restarting and 
APPENDing the eigenvector file. 
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9. 


USER MESSAGES AND OPTIONAL DIAGNOSTICS 


9.1 NASTRAN Functional Module User Messages for the Tridiagonal Reduction 
Method 

9.1.1 General 

The following is a description of the NASTRAN user messages which may be 
generated by NASTRAN during the execution of the Tridiagonal Reduction method 
and which are unique to this method. Explanatory information is provided follow- 
ing the text of each message and, in the case of a fatal message, corrective 
action is indicated. Refer to the NASTRAN Users’ Manual, Section 6 for a com- 
plete listing of other system and user messages. 

Fatal messages cause the termination of the execution following the print- 
ing of the message text. These messages will always appear at the end of the 
NASTRAN output. Warning and information messages will appear at various places 
in the output stream. Such messages convey only warnings or information to the 
user. Consequently, the execution continues in a normal manner following the 
printing of the message text. 

9.1.2 List of User Messages 

2385*** USER WARNING MESSAGE 2385, DESIRED NUMBER 0F EIGENVALUES EXCEED 
THE EXISTING NUMBER, ALL EIGENS0LUTI0NS WILL BE S0UGHT . 

The desired number of eigenvalues specified on the EIGB card 
(NEP) or the EIGR card (ND) exceeds the rank of the [KS| a ] or [M ] 
matrix, which is the maximum number of existing eigenvalues. 

2386*** USER FATAL MESSAGE 2386, STIFFNESS MATRIX SINGULARITY CANN0T BE 
REM0VED BY SHIFTING. 

Check the specification of masses on C0NM1, C0NM2, CMASSi, 
material defintion and element property cards to ensure that the 
degrees-of-freedom in the analysis set are not all massless. 

2387*** USER WARNING MESSAGE 2387, PR0BLEM SIZE REDUCED T0 **** DUE T0 
0RTH0G0NALITY DRIFT 0R NULL TRIAL VECTOR. 

ALL EXISTING MODES MAY HAVE BEEN 0BTAINED . USE DIAG 16 T0 
DETERMINE ERR0R B0UNDS . 
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2388*** 

2389*** 


2390*** 


2391*** 

2392*** 


The Tridiagonal Reduction Method cannot generate a reduced 
problem size of the order prescribed in Section 10,6.2.3 of 
the Theoretical Manual. However , the desired number of accurate 
eigenvalues specified on the EIGB card (NEP) or the EIGR card 
(ND) may have been obtained. A detailed list of the computed 
error bounds can be obtained by requesting DIAG 16 in the 
EXECUTIVE C0NTR0L DECK. 

USER WARNING MESSAGE 2388, USER SPECIFIED RANGE N0T USED F0R 
FEER BUCKLING, THE R00TS 0F L0WEST MANGITUDE ARE 0BTAINED. 

The value of LI specified on the EIGB card is ignored for 
buckling analysis by the Tridiagonal Reduction (FEER) method. 

USER WARNING MESSAGE 2389, PR0BLEM SIZE REDUCED. N0 M0RE TRIAL 
VECT0RS CAN BE 0BTAINED . 

The desired number of eigenvalues specified on the EIGB card 
(NEP) or the EIGR card (ND) exceeds the number that can be 
calculated by the Tridiagonal Reduction (FEER) method. Check 
whether the requested number of eigenvalues exceeds the rank 
of the tK| a ] or [M aa ] matrix, which equals the number of exist- 
ing eigenvalues. 

USER WARNING MESSAGE 2390, **** FEWER ACCURATE EIGENS0LUTI0NS 
THAN THE **** REQUESTED HAVE BEEN FOUND. USE DIAG 16 T0 
DETERMINE ERR0R B0UNDS . 

The number of eigenvalues passing the eigenvalue relative-error 
test is less than the number requested on the EIGB or EIGR card. 
The maximum allowable error is specified in field 5 on the above 
cards. A detailed list of the computed error bounds can be ob- 
tained by requesting DIAG 16 in the EXECUTIVE C0NTROL DECK. A 
checkpoint and restart should be employed to obtain additional 
accurate eigensolutions. 

USER FATAL MESSAGE 2391, PR0GRAM L0GIC ERROR IN FEER. 

USER INF0RMATI0N MESSAGE 2392, **** M0RE ACCURATE EIGENS0LUTI0NS 
THAN THE **** REQUESTED HAVE BEEN FOUND. USE DIAG 16 T0 DETERMINE 
ERR0R B0UNDS. 

The number of eigenvalues passing the eigenvalue relative-error 
test is greater than the number requested on the EIGB or EIGR 
card. The maximum allowable error is specified in field 5 on 
the above cards. A detailed list of the computed error bounds 
can be obtained by requesting DIAG 16 in the EXECUTIVE C0NTR0L 
DECK. 
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2393*** 


USER WARNING MESSAGE 2393, THE REDUCED- SYSTEM EIGENVECTOR 
CORRESPONDING T0 EIGENVALUE **** D0ES N0T MEET CONVERGENCE 
CRITERION. ABS0LUTE RELATIVE ERROR BETWEEN SUCCESSIVE 
ITERATES IS ****.. 

The accuracy of the corresponding physical eigenvector is in 
doubt. Refer to the Eigenvalue Summary Table for the largest 
error in the generalized mass matrix. 

9.2 The Eigenvalue Summary Table 

The following summary of the eigenvalue analysis performed is automatically 
printed when rigid formats using the Tridiagonal Reduction (FEER) method are in- 
voked : 

1. Number of eigenvalues extracted. 

2. Number of starting points used. 

This corresponds to the total number of random starting and re- 
start vectors used by the FEER process. 

3. Number of starting point moves. 

Not used in FEER (set equal to zero). 

4. Number of triangular decompositions. 

Always equal to one, except for unshifted vibration problems 
(roots starting from the lowest requested) . In this case a 
maximum of three shifts and three decompositions are employed 
to remove possible stiffness matrix singularities. 

5. Total number of vector iterations. 

The. total number of reorthogonalizations of all the trial 
vectors employed. 

6. Reason for termination. 

(0) Normal termination 

(1) Fewer than the requested number of eigenvalues and eigen- 
vectors have been extracted. 

(3) The problem size has been reduced. However, the desired 
number of accurate eigensolutions specified on the EIGB 
or EIGR card may have been obtained. A detailed list of 
the computed error bounds can be obtained by requesting 
DIAG 16 in the EXECUTIVE C0NTROL DECK. 
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7. Largest off-diagonal modal mass term and the number failing the 
mass orthogonality criterion. 

9.3 Optional Diagnostic Output 

The user can obtain special detailed information relating to the genera- 
tion of the reduced problem size, the elements of the reduced tridiagonal matrix, 
computed error bounds and other numerical tests by requesting DIAG 16 in the 
NASTRAN Executive Control Deck. 

The meaning of this information is explained below in the order in which 
it appears in the DIAG 16 output. 

0RDER - The order of the unreduced problem (size of the [K ] matrix) 

aa 

MAX RANK - The maximum number of existing finite eigensolutions as initially 

detected by FEER 

RED 0RDER - The order of the reduced eigenproblem which will be solved to 
obtain the number of accurate solutions requested by the user 

0RTK VCT — The number of previously computed accurate eigenvectors on the 

eigenvector file which were generated prior to a restart or by 
the NASTRAN rigid body mode generator 

USER SHIFT - Used only in frequency problems. The user specified shift after 
conversion from cycles to radians - squared 

INTERNAL SHIFT- Used only in frequency problems. A small positive value auto- 
matically computed to remove singularities if the user has 
specified a zero shift. Otherwise, the negative of the user 
shift 

SINGULARITY CHECK - PASS: the shifted stiffness matrix is non-singular 

****: the number of internal shifts needed to remove 

stiffness matrix singularities 

TRIDIAG0NAL ELEMENTS R0W j, **, ***, **** - Lists the computed tridiagonal ele- 
ments of the reduced eigenmatrix: 

j - Matrix row 
** - Diagonal element 

*** - Off-diagonal element 

**** - First estimate of off-diagonal element in the next row 
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0RTH ITER 


The number of times a reorthogonalization of a trial vector 
has been performed. 


MAX PR0J 


N0RMAL FACT 


The maximum projection of the above trial vector on the pre- 
viously computed accurate trial vectors (prior to the current 
reorthogonalization) 

The normalization factor for the reorthogonalized trial vector. 


0PEN C0RE N0T USED *** FEER3 - open core not used by Subroutine FEER3, in 
single-precision words 

FEER QRW ELEMENT *, ITER **, ***, RATIO ****, PR0J *****; 

* - The internal eigenvalue number in the order of its extraction by 

FEER 


** - The number of inverse power iterations performed to extract the 

associated eigenvector of the reduced system (this is not a 
physical eigenvector) 

*** - If a multiple root has been detected, the number of times that 

the previous multiple-root, reduced-system eigenvectors have been 
projected out of the current multiple-root eigenvector before re- 
peating the inverse power iterations 

**** - The absolute ratio of maximum, reduced-system eigenvector elements 

for successive inverse power iterations 

***** _ xhe maximum projection of a current multiple-root eigenvector on 

previously computed eigenvectors for the same root. 


PHYSICAL EIGENVALUE *, **, THE0R ERR0R *** PERCENT, PASS OR FAIL: 


* 


The internal eigenvalue number in the order of its extraction by 
FEER 


** - The associated physical eigenvalue (X for buckling problems, 

0)2 for frequency problems) 

*** - Theoretical upper bound on the relative eigenvalue error, in 

percent 

PASS - The computed error is less than or equal to the allowable specified 

on the EIGB or EIGR bulk data card (default is .001% where n is the 
order of the stiffness matrix) 

FAIL 

The computed error is greater than the allowable and this mode is 
not accepted for further processing 


0 PEN C0RE N0T USED *** FEER4 - open core not used by Subroutine FEER4 , in single 
precision words 
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FEER COMPLETE *, **, ***, **** 


* 


** 

*** 


The remaining CPU time available following decomposition of the shifted 
stiffness matrix, in seconds (the total time is specified on the TIME 
card in the Executive Control Deck) 

The remaining CPU time, in seconds after completing Subroutine FEER3 

The remaining CPU time, in seconds after completing Subroutine FEER4 

The total operation count for FEER after decomposition of the shifted 
stiffness matrix. One operation is considered to be a multiplication 
or division followed by an addition 
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